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Since the 1980s it has been possible to probe crystalhzed matter, thanks to X-ray 
or neutron scattering techniques, to obtain an accurate charge density or spin 
distribution at the atomic scale. Despite the description of the same physical 
quantity (electron density) and tremendous development of sources, detectors, 
data treatment software etc., these different techniques evolved separately with 
one model per experiment. However, a breakthrough was recently made by the 
development of a common model in order to combine information coming from 
all these different experiments. Here we report the first experimental 
determination of spin-resolved electron density obtained by a combined 
treatment of X-ray, neutron and polarized neutron diffraction data. These 
experimental spin up and spin down densities compare very well with density 
functional theory (DFT) calculations and also confirm a theoretical prediction 
made in 1985 which claims that majority spin electrons should have a more 
contracted distribution around the nucleus than minority spin electrons. 
Topological analysis of the resulting experimental spin-resolved electron density 
is also briefly discussed. 



1. Introduction 
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We learn from quantum physics that, in a crystal or a molecule, 
electrons are smoothly distributed according to a probability 
distribution law. In a bonded system, this electron distribution 
is significantly different from that of the mere sum of non- 
interacting atoms (the 'independent atoms model', lAM). This 
difference is indeed the key quantity to understand (and 
predict) not only chemical bonding mechanisms but also the 
resulting properties of molecules and crystals. When the 
system under study possesses magnetic properties, the charge 
distribution alone is not sufficient. One needs then to gain 
better knowledge of how the electrons, depending on their 
spin state, spread around and between the nuclei. 

The purpose of this article is to demonstrate that the 
experimental determination of the electron position prob- 
ability density has now reached a point where it has become 
possible to model and separate up and down spin carrier 
contributions to the distribution. However, this type of result 
is only possible by combining different types of experiments 
and interpreting them through a common model (Deutsch et 
al, 2012). In this paper we report on the first successful 
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experimental reconstruction of a spin-resolved electron 
density distribution and some insights obtained in this way for 
a model molecular magnetic crystal. 

2. X-rays for an electron probability density distribution 

Since the pioneering experiments conducted by Laue a 
hundred years ago it has been established that an X-ray 
diffraction (XRD) experiment provides a measurement of the 
Fourier coefficients (named 'structure factors', ^xrd(Q)) 
the electron distribution in the crystal. While the Fourier 
decomposition is infinite, the accuracy of an electron density 
reconstruction is limited by the experimental resolution ^in 
reciprocal space, typically around Q^^^ = sin^/A = l.lA . 

Owing to such resolution limitations and the need for 
experimental corrections, the electron distribution is often 
reconstructed by means of a model. The most popular one is 
the 'deformed pseudo-atoms' model proposed by Hansen & 
Coppens (1978). This model assumes that the electron prob- 
ability density distribution, denoted hereafter as p(r), can be 
fairly expressed as a sum of densities centred on each nucleus 
- hence the name 'pseudo-atoms'. Initially, each of these 
atomic contributions is identical to what it would be if the 
atoms were isolated. In particular, the atomic electron density 
is assumed to be locally isotropic and neutral. This is indeed a 
common hypothesis that is always used for crystal structure 
determination. However, upon chemical bond formation and 
subsequent charge transfers, the outer part of these pseudo- 
atomic electron densities wiU need to expand (or contract) and 
distort in some specific directions. For each of these pseudo- 
atoms, an expansion over a limited set of real spherical 
harmonics functions y£^±(^, (p) has been found to efficiently 
account for the amount of information provided by X-ray 
diffraction experiments. Therefore, the global contribution of 
a given atom / to the density distribution within the Hansen- 
Coppens (H-C) model is 

p<«'-)(r)+p(™'>^pr(^,r) 

I m=0 

where Rff\K'ir) are radial Slater- type functions. 

The best suited values for extension/contraction {/c^, a:-} and 
valence and multipole population {P^i^^^\ Pix,m±) parameters 
are determined to match experimental data {^xrd(Q)}- 

This particular model for charge density has encountered 
such tremendous success since the 1980s that high-resolution 
XRD is now frequently employed to analyze chemical bonds 
and more generally the charge redistribution due to molecule 
or crystal formation. This experimental approach to electron 
probabihty distribution has become extensively used and was 
developed to calibrate advanced quantum chemistry compu- 
tations. Such a landmark is of course essential, especially if the 
aim is to extend the theoretical modelling to a set of systems 
which cannot yet be studied experimentally (Gatti & Macchi, 
2012). 



3. Polarized neutrons for unpaired electron probability 
distribution 

As a consequence of this success, the quest for a deeper 
understanding of mechanisms responsible for spin magnetism 
has led to the use of the same approach to reconstruct electron 
magnetization distribution in crystals from polarized neutron 
diffraction (PND) data. 

While XRD and non-polarized neutron diffraction data 
consist of integrated intensities of Bragg reflections, PND 
measures 'flipping ratios' (hereafter denoted ^pnd(Q))- They 
are defined as the ratio between the diffracted intensities for 
spin up and spin down incident neutrons. This technique thus 
aims to identify regions of the crystal unit cell where an 
electron spin majority exists and how it is distributed. PND 
gives access to the magnetization density that is the sum of the 
pure spin density and the orbital contribution (Schweizer, 
2006). 

In a PND data collection, the explored reciprocal space 
domain is generally more limited than its XRD counterpart 
because of the restricted access to neutron beam time and the 
available range of experimental geometries (the sample is 
submitted to a strong magnetic field, typically 5 T). Thus, PND 
is often limited to the domain of Q for which a significant 
magnetic contribution is expected, i.e. Q^ax is 0.5 for 3d 
transition metals. Moreover, as opposed to structural data 
collections which concern all the existing reflections in a given 
domain of Q, PND measurements are performed only for the 
strong nuclear reflections for which good experimental 
accuracies may be achieved for the flipping ratios. Therefore, 
only ~ 10% of all reflections in a given reciprocal domain are 
measured by PND. 

Despite experimental difficulties that hinder its resolution, 
PND has provided a wealth of information to the magnetism 
community. In fact, the spin density visualization provides 
direct information about the magnetic interaction pathways 
and the sign of the spin density distribution along this 
pathway, i.e. a crucial test of the nature of the mechanism 
(Pontillon et al, 1999; Pillet et al, 2001). 

4. The combination challenge in determining a spin- 
resolved electron density distribution 

From the above considerations it is quite clear that XRD and 
PND consider electron distributions from different and 
complementary perspectives. XRD enables the reconstruction 
of the total electron distribution, p(r), while PND provides 
information which yields the spin density, ^(r). The two 
quantities can therefore be expressed in terms of spin-resolved 
electron densities, p^(r) and P|(r), respectively, representing 
spin up and spin down only electron density distributions 

p(r) = p^(r) + p^(r) (2) 

^(r) = p^(r) - p^(r). (3) 

It is therefore obvious that a combined analysis of accurate 
high-resolution X-ray and polarized neutron diffraction data 
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should yield unprecedented access to spin-resolved electron 
densities for crystals with significant magnetic properties. 
Despite such an elementary relationship between spin and 
charge densities no successful combined analysis has been 
reported so far. The reason mostly lies in the lack of a common 
model and the difficulty in finding a fair balance between sets 
of data obtained from different experiments. 

4.1 . A 'spin-split' pseudo-atoms model 

In order to carry out the present work, a 'spin-split' pseudo- 
atoms model was adopted to simultaneously interpret XRD, 
unpolarized neutron and PND sets of data (Deutsch et al, 
2012). Following the above mentioned H-C model, each 
atomic contribution is modeled by a frozen core density 
(identical to that of a similar free atom) and a parameterized 
valence allowing for contraction/expansion as well as angular 
distortion with a similar formalism, as explained above. 
However (and thereby departing from the usual H-C model), 
when an atomic region is expected to carry a significant spin 
distribution, the valence part of the model is split into spin up 
and spin down contributions. The number of parameters in the 
model is thus moderately increased compared with the usual 
charge density model for the total density distribution as the 
magnetization distribution is often limited to a small set of 
atoms. For an atom expected to carry a magnetic moment, the 
pseudo-atom contribution p(r) then writes 

I m=0 
I m=0 

where and refer to spin up and spin down parameters, 
respectively. Thus, the challenge is in the determination of P^ 
and P^, as well as /c^ and /c^, against XRD and PND data in a 
unique refinement procedure, thanks to the complementa- 
rities of the two methods. 

Finding a fair way to gather experimental data from 
different origins still represents a real challenge. No routine 
strategy has yet been identified and each specific case has to be 
carefully considered. In our approach it was found that a 
satisfying result was reached when the following quantity was 
considered for minimization (Gillet et al, 2001; Bell et al, 
1996): 



.log 



" ' AKU 

E 



+ ^PND log 



E 



or 



(5) 



where A^xrd ^^^1 A^pnd ^^e the number of data points in each 
set and a represents the estimated standard deviation for each 
data point. 

A more detailed description of the joint refinement meth- 
odology is given in Deutsch et al. (2012). 



5. First spin-resolved electron density reconstruction 
and comparison with ab initio quantum computations 

5.1. A molecular magnet as a test case 

Crystallized molecular magnets are typical examples of 
large systems with a limited number of magnetic atoms. Thus, 
the molecular complex Cu2L2(N3)2 [L = l,l,l-trifluoro-7- 
((dimethylamino)-4-methyl-5-aza-3-hepten-2-onato)] recently 
studied by PND ( Aronica et al. , 2007) was chosen as a test case 
for this joint investigation. This complex belongs to an 
extensively studied family of dicopper complexes in which the 
Cu^^ ions are coupled by two azido bridges (Ns"). This family 
is of particular interest in the field of molecular magnetism 
because the nature of the intramolecular magnetic coupling 
varies from strongly ferromagnetic to strongly anti- 
ferromagnetic coupling depending on the geometry of the 
bridging ligands. In fact, when two azido groups are symme- 
trically connected to the Cu atoms by only one N atom (End- 
On) a triplet ground state is always observed, while a singlet 
state is observed when these two azide groups relate the Cu 
atoms by the extremities (End-to-End). 

For a long time, the Cu— Cu magnetic interaction 
mechanism was controversial. It was finally established, partly 
thanks to PND (Aebersold et al, 1998), that the ferromagnetic 
coupling in End-On (EO) complexes is due to the quasi- 
orthogonality of the magnetic orbitals centred on the Cu 
atoms for a specified range of values of the Cu— N — Cu 
bridging angle, while the antiferromagnetic (AF) coupling in 
End-to-End (EE) complexes originates from the strong 
overlap between these orbitals. 

In the case of asymmetric bridges, a more subtle interaction 
mechanism has to be invoked. Indeed, some EO compounds 
have been shown to exhibit AF coupling and, conversely, 
ferromagnetic coupling in EE compounds have recently been 
obtained such as in the present End-to-End Cu2L2(N3)2 
complex. The crystal structure is formed by discrete neutral 
centrosymmetric five-coordinated Cu^^ dinuclear complexes 
[Cu- • -Cu = 5.068 (1) A; Aronica et al, 2007], which are weU 
isolated from each other in the crystal. The azido groups 
bridge the Cu ions in an asymmetric fashion [Cu— N3 = 
2.000 (1) A, Cu-N5 = 2.346 (1) A; Fig. 1]. 

5.2. Results 

Details of the three distinct data collections on single 
crystals of the molecular complex Cu2L2(N3)2 by X-ray, 
neutron and polarized neutron diffraction, respectively, are 
reported in Table 1, as well as the conditions and statistical 
agreement factors of the joint refinement. For a detailed 
description of the first XRD and PND data treatment for the 
copper compound see Deutsch et al. (2013) and Lecomte et al. 
(2011). 

The deformation charge density upon the molecule bond 
formation displayed in Fig. 2(a) was deduced by subtracting 
the density of independent atoms (computed) from the 
experimental electron distribution Ap(r) = p^(r) + P|(r) — 
PiamW reconstructed from PND and XRD data. Fig. 2(b) 
shows the experimental spin density s(r) = p^r) — p^(r) which 
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Table 1 

Details of the data collections and combined refinement (the total number of parameters in the joint 
refinement is equal to 887, which is smaller than the sum of parameters refined on the different data sets 
because some parameters are common to different experiments). 





X-ray diffraction 


Neutron diffraction 


Polarized neutrons 


Number of measured reflections 


15 731 


5049 


474 




Number of used reflections 


7208 (/ > 3cr, sin Oik < 1 A" 


^) 2303 (/ > 3cr) 


212 (Fn 


> 5 X 10"^^ cm) 


Resolution (sin 6IX)^^^ (A~^) 


1.13 


0.78 


0.5 




Temperature (K) 


10 


30 


2 




Magnetic field (T) 






5 




Wavelength (A) 


0.71073 


0.832 


0.84 




Statistical agreement factors (%) 


R^{F) = 1.52 


R^{F^) = 5.50 




= 8.96 








RUR) = 


0.98 


Number of parameters 


626 


297 


69 





was obtained taking into account a correction for the orbital 
contribution to the total magnetization density. This correc- 
tion was done by applying the usual dipole approximation to 
describe the orbital form factor (Squires, 1978) which is 
written as a sum of the tabulated radial integrals (JO) 



End-On 



<> < > 



N 
I 

N 



End-to-End 




(b) (c) 
Figure 1 

Di-azido copper complexes. Schematic representation of (a) End-On and 
(b) End-to-End conformation of di-azido di-Cu complexes, (c) View of 
the Cu2L2(N3)2 molecule. N atoms are represented in blue, O in red, C in 
grey, F in yellow and Cu in orange. H atoms are not shown for reasons of 
clarity. 




Figure 2 

Charge and spin density maps in the plane containing Cu, Ol and N5. (a) Static deformation density map 
obtained by means of the joint refinement strategy. Isocontours are drawn for dz 0.01 x 2" e A~^ with n = 
0-13 (positive red, negative blue), (b) Spin density map obtained by means of the joint refinement 
strategy. Isocontours are drawn for ± 0.01 x 2" /xB A~^ with n = 0-13, spin up contours in red, spin down 
contours in blue. 



and (jl) for Cu^^ (Brown, 1992) 
multiplied by a coefficient = 
^sfecu2)/gcu. where gcu = 2.175 
(Aronica et al, 2007). A value of 
0.07 |iB was taken for this coeffi- 
cient, according to the refined 
valence population associated with 
the orbital form factor in the PND- 
only refinement. 

The charge deformation density 
map shows the usual distribution of 
the electrons in all bonds. It 
describes the d electron redis- 
tribution around the Cu atom, and 
the lone-pair densities of N and O 
atoms of the ligand directed toward the Cu atom. The corre- 
sponding spin density exhibits a (i^2_^2-type distribution 
around the copper with maxima directed toward the N and O 
atoms and depletion in the bisecting direction. Very httle spin 
density is found in the vicinity of the ligand atoms. 

The first spin-resolved valence electron density distribu- 
tions obtained experimentally from a joint analysis of a 
common model are displayed in Fig. 3. As can be clearly 
observed, the spin-split model obtained from the proposed 
joint refinement strategy has made it possible to successfully 
discriminate the density probability distribution of spin up and 
spin down electrons. Quantum computations (Frisch et al., 
2009), by means of the Density Functional Theory (B3LYP/ 6- 
31-h+Gdp), were carried out on an isolated molecule in its 
experimental geometry. The theoretical and experimental 
distributions compare extremely well (Fig. 3); the differences 
are hardly seen at the drawn contour level. The spin up 
distribution in the vicinity of the copper nucleus is spherical, 
while the down spin distribution shows maxima in the 
bisecting direction of the ligands in both cases (theory and 
experiments). It thus appears 
striking that most of the electron 
anisotropy around the Cu atom 
(Fig. 3) should be attributed to spin 
down electrons. This is indeed 
confirmed by a J-type function 
analysis (Holladay et al., 1983) 
which neglects the covalent part of 
the metal-ligand bonds of the Cu 
atom reported in Fig. 4. It is found 
that 30% of spin down electrons he 
in the J^^-type function with 
corresponding d^2_y2 depletion 
(9%), while all d^^, dy^ and d^i are 
almost equally populated. 

One utmost consequence of the 
spin-resolved model is that it is 
shown for the first time that the 
valence spin f density is 5% more 
contracted than the spin f density 
[/ct = 0.998 (1), = 0.943 (1)]. 
This is in agreement with theore- 
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Figure 3 

Spin-resolved electron densities. Left: (a) Experimental spin up (majority) and (c) experimental spin 
down (minority) valence electron densities from joint refinement of the spin-split model. Right: (b) 
Theoretical spin up (majority) and (d) theoretical spin down (minority) valence electron densities from 
ab initio quantum computation. The density distributions are represented in the Cu— Nl — Ol plane 
(contours 0.01 x 2" e A"^ (n = 0-12)). 



tical predictions (Watson & Freeman, 1960) based on a /c 
refinement of calculated spin-resolved magnetic form factors 
(Becker & Coppens, 1985) which indicated a contraction of 
1% of the electronic cloud with spin f . 

The topological analysis of the experimental electron 
density provides a precise description of the interactions in the 
crystal (Bader, 1990). It locates particular points, called critical 
points (CPs), where the gradient of the density vanishes. There 
are four types of critical points: local maximum of the density. 



30% 



19% 21% 20% 18% 21% 23% 20% 20% 

li 11 li li 



xz 



yz 



xy 



Figure 4 

Schematic representation of the Cu J-orbital type function populations 
for up and down electrons. The arrows sizes are proportional to the 
respective spin populations. 



local minimum of the density and 
two saddle points. The saddle 
points with two negative curvatures 
are called bond critical points 
(BCPs), among them two kinds can 
be distinguished by the Laplacian 
value [V^p(r)]: 'shared-shell' inter- 
action, where p(CP) is high and 
V^p(CP) < 0 which means a 
concentration of electron density at 
the BCP, and 'closed-shell' inter- 
action, where p(CP) is low and V^p 
(CP) > 0 which means a depletion 
of the electron density at the BCP. 
The present work using XRD and 
PND joint refinement allows for 
the first time the computation of 
experimental spin-resolved topolo- 
gical properties (Souhassou & 
Blessing, 1999). To the best of our 
knowledge there is only one paper 
devoted to spin-resolved topolo- 
gical analysis which is a theoretical 
analysis of b.c.c. and f.c.c. iron 
(Jones et al., 2008). The features of 
the critical points around the Cu 
atoms are summarized in Table 2 
and the spin-resolved Laplacian 
maps are displayed in Fig. 5. The 
BCPs for the Cu environment are 
of 'closed-shell' type with a small 
value of the electron density at the 
BCP. The large positive value of the 
Laplacian indicates dominant ionic 
interactions. The positions of the 
BCPs are almost equal whatever 
the spin: p(CP) and V^p(CP) are slightly higher for spin up 
electrons, which is both due to a more contracted spin up 
density and a higher number of spin up electrons (0.74 e 
more); however, considering the uncertainty on these values 
(around 10% for the Laplacian) the differences may not be 
significant. 

Moreover, a good agreement between theory and experi- 
ment is observed for the net charges as obtained by integrating 
spin up or down electrons over the corresponding atomic 
basins for the Cu atom (Table 3). The majority valence spin 
state (t) is experimentally 0.74 e greater than the minority 
one (I) compared with the DFT calculation, 0.63 e. 

6. Conclusion 

It clearly appears that the 'spin split' model, together with the 
proposed joint refinement method, is successful at precisely 
retrieving all the essential features of the electron distribution 
with their respective spin dependences. Most important is the 
dramatic difference between spin up and spin down angular 
distributions in the vicinity of the copper nuclei. The spin- 
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Table 2 

Topological properties of the spin up and spin down densities at critical 
points in the neighbouring Cu atom. 



Table 3 

Experimental copper spin up (^f) and spin down (q^) electron 
populations (in electrons) integrated over Bader atomic basins. 



Bond 


d(A) 


(a A ~3\ 

(e A ) 


V pKyv) 
(e A ) 


"l-bcp 

( °S.\ 
(A) 


1 

(e A ) 


1 

(e A ) 


1 

(e A ) 


Spin up 
















Cu-Ol 


1.936 


0.36 


6.38 


0.966 


-2.05 


-1.90 


10.33 


Cu-Nl 


1.961 


0.34 


5.90 


0.979 


-1.65 


-1.53 


9.09 


Cu-N5 


2.003 


0.31 


4.51 


1.002 


-1.57 


-1.52 


7.60 


Cu-N2 


2.036 


0.30 


4.56 


1.004 


-1.57 


-1.34 


7.47 


Spin down 














Cu-Ol 


1.936 


0.33 


5.73 


0.965 


-1.83 


-1.66 


9.22 


Cu-Nl 


1.960 


0.33 


5.48 


0.970 


-1.59 


-1.43 


8.50 


Cu-N5 


2.003 


0.29 


4.09 


0.997 


-1.47 


-1.39 


6.94 


Cu-N2 


2.037 


0.29 


4.15 


1.000 


-1.47 


-1.21 


6.84 



resolved electron density model confirms for the first time 
experimentally the different contraction of spin up and spin 
down radial distribution of Cu atoms as predicted by theory. 
The topological analysis of the two distributions shows a very 
good agreement between theory and experiments. We have 
demonstrated that the joint refinement gives access to the spin 
up and spin down distributions that were not accessible by 
standard refinements. This can be applied to any magnetic 
crystalline material, inorganic or organic. In a wider 
perspective, this method paves the way for combining 
different scattering or spectroscopic experiments as electron 
diffraction, CBED, NQR, NMR, EPR . . . , that will lead to a 
more precise and robust description of the electronic behavior 
in crystalline solids. 
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q\ (e) 


qi (e) 


(e) 


Experimental 


14.55 


13.81 


0.74 


Theoretical 


14.29 


13.66 


0.63 




Figure 5 

Spin-resolved Laplacian maps in the plane containing Cu, Ol and N5: ia) spin up, {h) spin down; b.c. for 
bond critical point and r.c. for ring critical point (saddle point with two positive curvatures). 
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